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1 Introduction 

This article deals with applications of the Moutard transformation |T2] which 
is a two-dimensional version of the well-known in solitonics Darboux trans- 
formation to some problems of the spectral theory of two-dimensional oper- 
ators and of (2 -|- l)-dimensional nonlinear evolution equations. 

In particular, the main results consist in the explicit construction of 

• two-dimensional Schrddinger operators 

H = -A + u = -{dl + dl) + u{x, y) 

with fast decaying smooth rational potentials such that their L2-ker- 
nels are nontrivial and, moreover, contain at least two-dimensional 
subspaces spanned by rational eigenfunctions (see Theorems^and^; 

• blow-up of solutions to the Novikov-Veselov (NV) equation, which 
is a two-dimensional generalization of the Korteweg-de Vries (KdV) 
equation, with fast decaying rational Cauchy data (see Theorem^. 

The first construction was already announced and briefly sketched in |19) . 
For operators with such fast decaying potentials there exists a nice spectral 
theory [H [16]. We note that for one-dimensional Schrodinger operators 
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with fast decaying potentials the existence of square-summable eigenfunc- 
tions at zero energy level is impossible (see, for instance [S]) and for higher- 
dimensional operators (i.e., for dimensions greater or equal than 5) one may 
easily construct such examples for which the kernel contains a smoothing of 
the Green function of A. However for two-dimensional operators these are 
the first known examples with nontrivial L2-kernel. 
The Novikov-Veselov equation has the form 

Ut = d^U + B^U + 3d{VU) + 3d{VU) =0, 

Bv = dU, 

and it is the first in the hierarchy of equations which have the form 

Ht = HA + BH, 

where A and B are differential operators. Here and everywhere below we 
use the standard derivatives d = dz = ^{Bx — iBy), B = Bz = ^{Bx + iBy) 
w.r.t. z = X + iy. These equations were introduced in J7] as equations 
which preserve the zero-level spectrum of H. The principal part of the n-th 
equation from the NV hierarchy takes the form 

where . . . stays for terms of lower order. 

It follows from the inverse scattering method that solutions to the KdV 
equation with analytical fast decaying Cauchy data do not blow up and 
the well-posedness of the Cauchy problem for this equation is established 
in many functional spaces (see [13] and references therein). However as we 
show this is not true for its natural two-dimensional generalization. 

The examples of blow-up solutions were obtained not by using the inverse 
scattering method which is not well-developed in this situation. Indeed, the 
inverse spectral problem for a two-dimensional Schrodinger operator at a 
fixed energy level was first posed in |6j and has been studied for positive 
energy levels |T0] or levels below the ground state [H]. We think that the 
study of this problem on the zero energy level will help to understand both 
phenomena which we discuss in this article. 

We also have to remark that our potentials obtained by iterations of the 
Moutard transformation may be considered as two-dimensional generaliza- 
tions of one-dimensional rational solitons obtained in the same way using 



the Darboux transformation (see { 2.2). However in the one-dimensional case 



these potentials are always singular. 
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Finally we remark that our potentials are constructed by iterations of the 
Moutard transformation from the constant potential and all such potentials 
are integrable on the zero energy level in the sense that all solutions to the 
equation Htp = may be explicitly constructed via quadratures from linear 
combinations of harmonic functions. We give the details of this construction 
in Section [H 

The authors thank P.G. Grinevich and S.P. Novikov for useful discus- 
sions. 

2 Darboux and Moutard transformations 

2.1 The Darboux transformation 

Let 

be a one-dimensional Schrodinger operator and let a; satisfy the equation 

Huj = 0. 

The function uo determines a factorization of H: 

H = A^A, A = -^ + v, A'^ = ^+v, v = '^. (2) 
ax ax UJ 

Indeed we have 

A^A =(^ + ,\(-^+i\=-^+v^ + v^ 
\ax J \ ax J dx^ 

and the equation 

is equivalent to Hu) = 0. If w is real- valued we have A* = AJ . 

The Darboux transformation of H [5 is the swapping of AJ and A: 

H = A^ A^ H = AA'^ , 

or in terms of u: 

u = + Vx ^ u = — Vx- 
It is easy to check the following: 

Proposition 1 If (p satisfies the equation Hip = Eip with E = const then 
if = Aip satisfies the equation Hip = Eip. 

Remark 1. In general the Darboux transformation is defined for any 
solution to the equation Huj = civ with c = const. In this case it reduces to 
the transformation of H' = H — c for which H'lo = 0. 
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2.2 One-dimensional solitons via the Darboux transforma- 
tion 

Let u = and lo = uJi = x = Ti. Then 

1 1 ^2 

V=-, Vx = Ul = U=^. 

The function 

is meromorphic m. P = {E,\) on the Riemann surface T = {A^ = and 
for every E its branches give a basis of solutions to the equation 

Hii^ = (^-^+ui^%b = Ei; 

(we normahze ip by the condition ^/^ ~ ^i\/Ex as — > c>o). Now we may apply 
the Darboux transformation defined by a;2 = a;^ + ^ and obtain another 
potential. In particular, for every n the potential 

_ n(n + 1) 

is obtained after n iterations. The orbit of Un under the Korteweg-de Vries 
hierarchy is the n-dimensional family M„ of potentials. 
We have [DE]: 

1. there is a general recursion procedure found by Adler and Moser [1] 
for deriving the polynomials 

6'„(ri, . . . ,r„) 

such that 

(a) dn is a polynomial of degree "^'^2^"^'' in x = ri. In particular, we 
have 

Q\=X = Ti, 

(92 = + ra, 
6»3 = x^ + 5t2X^ + rax - 5r|; 
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(b) the function 

Un{x) = -2-^log6'„(ri + x,r2, . . . ,r„) 

is the potential of H obtained after n iterations of the Darboux 
transformation (which starts at uq = 0). Here Tk,k = 2, . . . , are 
free scalar parameters (integration constants) which appear one 
by one at every step of the iteration; 

(c) for ifn = %^ we have 

+ n„ = 0; 



2. Mn is an n-dimensional family parameterized by r G C" and consists 
of the potentials M„(ri + r2, . . . , r,i); 

3. there are birational transformations r ^ i of the form 

tk = cLkTk + gkin, ■ ■ ■ ,Tk-i), a/c / 0, 

such that 

Wk = ^^'(""^ 

is the k-ih KdV flow; 

4. the zeroes xi, . . . ,a;„(„+i)/2 of the polynomials On are evolved by the 
KdV-flows as some integrable Hamiltonian systems and, for instance, 
for the original KdV equation 

their dynamics is described by the Calogero-Moser system. 
2.3 The Moutard transformation 

Let be a two-dimensional potential Schrodinger operator and let a; be a 
solution to the equation 

Hlj = (-A + u)uj = 0. 
Then the (elliptic) Moutard transformation of H (cf. [15 J is defined as 

H = -A + u- 2Aloguj = -A-u + 2 ^ ^ . 
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Proposition 2 // 99 satisfies the equation Hip = 0, then the function 9 
defined from the system 



M), = -u;2/'n ^ (^^e)y=u;^'^) (3) 

\UJ/y \iVJx 



satisfies H6 = 0. 

Obviously, if 9 satisfies Q tlien 



C 

9+—, C = const, (4) 

satisfies ([S]) for any constant C. 

We shall use the following notation for the Moutard transformation: 

M^{u) = u = u-2Aloguj, M^{ip) = {9 + -, CgC}. 

In the one-dimensional limit the Moutard transformation reduces to the 
Darboux transformation. Indeed, let u = u{x) depend on x only and to = 
f{x)e^y. Then / satisfies the one-dimensional Schrodinger equation 

and the Moutard transformation reduces to the Darboux transformation of 
Hq defined by /: 

H — Hq — —7 — > H — Ho — —7 . 

li g = g{x) satisfies Hog = Eg, then Hip = with ip = e^yg{x). We derive 
from (3) that 9 = e^yh[x) satisfies H9 = d ii h = (^^ - ■^^ 5, 

i.e. /i is a multiple of the Darboux transform of h = — ^^^ Ag where 
Hq — c = A is the factorization of Ho — c defined by /. The inverse 
Darboux transformation is given hy g = ^ ^ ^ ((S ^ ^• 

Remark 2 . There is another two-dimensional generalization of the Dar- 
boux transformation called the Laplace transformation. It is defined for a 
more general operator, i.e. for the Schrodinger operator with an electromag- 
netic field: 

i/ = 4(a + /3)(-a + a)+u. 
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and has the form 

H = 4u{-d + a)u-^{d + p)+u. 

So if Hip = 0, then = {—d + a)ip satisfies the equation Hip = 0. In the 
one-dimensional hmit it also reduces to the Darboux transformation. Its 
relation to integrable systems was recently studied in |18] . 

3 Soliton potentials via the Moutard transforma- 
tion 

In [T9j we gave simple examples of fast decaying rational potentials of the 
two-dimensional Schrodinger operator which have a degenerate L2-kernel. 
These examples are constructed using the Moutard transformation as fol- 
lows. 

Main construction. Let 

Ho = -A = -A + uq 
be an operator with a potential uq = and let uji and uj2 satisfy the equations 

HqUJi = HqUJ2 = 0. 

We take the Moutard transformations Mt^^ and Mo^j defined by uji and uj2 
and obtain the operators 

Hi = -A + ui, H2 = -A + U2 

where ui = M^^{uq),U2 = M^^{uq). By the construction, we have 

HiM^^{uJ2) = 0, H2M^^{oJi) = 0. 

Let us choose some function 

01 G M^,{lo2) 

and put 

62 = - — 01 E M^,{u{). 

L02 

These functions define the Moutard transformations of Hi and H2 and we 
obtain the operators H12 and H21 with the potentials 

ui2 = Me^{ui), U21 = Me^{u2). 

The following classical lemma is checked by a straightforward computation 
which we omit. 
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Lemma 1 1. ui2 = U21 = u, i.e. the diagram 

uo — > Ui 

U2 > U12 = U21 

where 61 G M^^^{ijJ2), ^2 = ^ Mx;2('^i)' ^-^ commutative; 

2. For V'l = ^ cf^f^ ^^2 = -ff^ we have 

Htpi = = 

where H = —A + u. 

We note that in this construction we have a free scalar parameter C (see 
Q) for the choice of 9i G M(^j(a;2)- This parameter can be used in some 
cases to build a non-singular potential u and functions ipi and '02 • 

Now let us apply this construction to the situation when 

Uo = 0. 

The desired formula for the potential u is given by the following theorem 
which is derived by simple and straightforward computations which we omit. 

Theorem 1 Let 



LVi = pi{z) + pi{z) , UJ2=P2{z) +P2{z), 

where pi andp2 are holomorphic functions of z. Let us consider the Moutard 
transformation of the operator Hq = —A defined by toi. Then the corre- 
sponding transformation ^ of UJ2 has the form 

Oi = ( {pip2 - P2P1) + I {{PiP2 - Pip'2)dz + {pip'2 - p'iP2)dz) ) . 

Therefore the second iteration defined by 61 gives us the operator H = —A+u 
with 

u = -2Aloga;i - 2A log 6*1 = -2A log(cJi6li) = 
- 2A log i (^{piP2 - P2P1) + j {{PiP2 - Pip'2)dz + {P1P2 - PiP2)dz)^ . (5) 
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The free scalar parameter which we mentioned above appears as the 
integration constant in ([5]). 

Let us apply this theorem to obtain some interesting examples of oper- 
ators. We consider hereafter only the cases when loi and uj2 are real- valued 
harmonic polynomials. 

Example 1. Let 

uJi=x + 2{x'^ -y^)+xy, 002 = x + y + ^(x^ - y^) + 5xy, 
which in terms of complex polynomials is written as 

^^2 ^ I/O t;-\2 l~^ 



p,{z)=^l--jz' + -, p^{z) = -{3-5i)z' + —z. (6) 

For some appropriate constant C in 61 we obtain 

5120(1 + 8x + 2y + ITx^ + 17y^) 



u 



[im + 4x2 + 4y2 + 16x3 + 4x2y + 16xy^ + 4y3 + I7(x2 + y2)2)2 

(7) 

5120|l-F(4-i)z|2 



(160-h|z|2|2-h(4-i)z|2)2 

and 



X -I- 2x2 + xy — 2y2 



160 + 4x2 _^ 4y2 _^ 16x3 + 4x2y + 16xy2 4y3 + 17(x2 y2)2 ■ 

1p2 = 



2x + 2y + 3x2 _^ ^Q^y _ 3y2 



(8) 



160 + 4x2 4y2 + 16x3 + Ax'^y + 16xy2 + 4y3 _^ 17(3.2 _^ ^2)2 
(here we simplify the expressions for ipi and 1^2 by multiplying them by some 



constant). On page 21 we put the graphs of this potential and the solutions 

i'l, ip2- 

Theorem 2 ([19j) The potential u given by ^ is smooth, rational, and 
decays like for r ^ 00. 

The functions ipi and ^/'2 given by ^ are smooth, rational, decay like 
l/r2 for r ^ 00 and span a two-dimensional space in the kernel of the 
operator L = -A -h n : L2(K^) ^ ^2(^2). 

Here and in the next example we put r = Y^x2~+~y2_ 
Example 2. Let us take for uji and L02 two harmonic polynomials of the 
third order: 



uji=x + ^ — + 2(-x3 - 3x^y + 3xy^ + y^), 
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— xy 2 

W2 = X + y H - 4(3x y -y ) 

1 5 

which in terms of complex polynomials take the form 

p,(.) = (,-l).3+(l + |i).'+i, p,W=2,.3+(l + A) (9) 

Then the potential u and the functions V'l and i/'2 take the form 
Fo{x,y) Fi{x,y) F2{x,y) 

*''w:y)- *''w:y) 

with 

Fo(x,y) = -1280000(25+20x-287x2+60x3+1800x^-30y-600xy-300x2y+ 

313y2 + 60xy2 + 3600x^^2 - 300y^ + ISOO/), 

G(x, y) = 40000 + lOOx^ + 40x^ - 387x^ + 40x^ + SOOx^ - GOx^y- 

SOOx^y - 200x^y + lOOy^ + 40xy2 + 26x2y2 + SOx^y^ + 2400x'^y2 - 60y^- 

800xy^ - 400x2y^ + 413y^ + 40xy*^ + 2400x2y^ - 200y^ + 800y^ 

y) = -lOx - 2x2 + 20x3 + 6xy + GOx^y + 2y2 - GOxy^ - 20y3, 

F2{x, y) = -lOx - 5x2 _ -^Qy ^ 2xy + 120x2y + 5y^ - 40y^ 

(as in the previous example we simplify the expressions for ipi and 1^2 mul- 
tiplying them by some constant). On pages 21 -22 we put the graphs of this 
potential u and the solutions ipi, ip2- 



Theorem 3 ([19j) The potential u given by (10) is smooth, rational, and 
decays like for r ^ 00. 

The functions "01 cLnd ip2 given by (10) are smooth, rational, decay like 
l/r^ for r ^ 00 and span a two-dimensional space in the kernel of the 
operator L = -A + n : ^2(^2) ^ L2(M2). 

Remark 3. We guess that for every > by applying this construction 
to other harmonic polynomials one can construct smooth rational potentials 
u and the eigenfunctions V'l and ip2 decaying faster than 
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4 Soliton equations 

4.1 Explicit solutions to the Novikov— Veselov equation 

For simplicity, in this section we renormalize the Schrodinger operator as 
follows 

1 V 

H = dd + U = -A-- (11) 
4 4 

with the standard d = dz = ^{d^ — idy), d = dg = ^{dx + idy). 

The Moutard transformation defined by a function u) takes the form 

U + 2dd\oguj 

and the transformation of eigenfunctions is given by the same formulas ([3| 
which are rewritten in terms of a complex coordinate as 

or 

(p ^ 9 = — / i^pdoj — LL!dip)dz — ((pdu> — ujdip)dz. 
^ J 

Let us assume that the function ip depends also on the time t and satisfies 
the equations 



Hip = 0, 

dtp = {d^ + d^ + 3Vd H 

where 



, -o (12) 
dtp = {d^ + d^ + 3Vd + 3V*d)ip 



dV = dU, dV* = dU. 
Then by straightforward computations we derive 

Proposition 3 (The extended Moutard transformation) The sys- 
tem (12) is invariant under the extended Moutard transformation 

ip —I- 9 = — / {pduj — ujdip)dz — {pBuj — ujdip)dz+ 



[pd^uj — ujd^p + ujB^p — pd^oj + 2{d'^pduj — dpd'^uj)— (^i^^^ 
2{d'^(pdui — dpB^uj) + 3V{pduj — ujdp) + 3V*{ijjdp — pduj)]dt, 
U + 2ddlogLu, V + 2d^logio, V* ^ V* + 28"^ log u . 



^ This theorem was first formulated in [T^. However in [14] the third fine in the formula 
for the Moutard transformation <p ^ 8 was omitted and here we make this correction. 
Due to this fact the formula (6.2.5) in [T3] does not give a solution of the Novikov- Veselov 
equation. Recently we learned that such a correction was already done in [12]. See also 
[3] for another approach to the extended Moutard transformations. 
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If uj is a real-valued function, then the latter transformations preserve the 
property V* = V. 



The compatibility condition for the system ( |12[ ) is the system 
Ut = d^U + d^U + ?,d{VU) + ?,d{V*U) = 0, 
dv = dU, dV* = Bu 

which for 

V* = v 

reduces to the Novikov-Veselov equation. In the sequel, we consider only 
the cases when the above equality holds. 

Let us consider the flow on the space of functions p{z, t) holomorphic in 

z: 

dp^d^ 

dt dz' ^ ^ 



If pi{z,t) and p2{z,t) satisfy (14) then the functions uj = pi + pi and 
If = P2 +p2 satisfy the system ( 12 ) with U = 1/ = 0. The extended Moutard 
transformation differs from the original one by the dt-term: 

^i{uj, ip)dz + ^'2(w, ip)dz + Q{u>, ip)dt 



Let us substitute uj = pi+pi, ip = P2+P2 into ^'5. By ( 14 ) and V = V* = 
the integrand is a closed form which implies 

+Pi,P2 +P2) _ dQ{pi +pi,P2 +P2) 
dt dz ' 

3^'2(pi +Pl,P2 +P2) dQ{pi+pi,P2+P2) 



dt dz 
Hence Qz and Qz are defined by and we have the following analog of 
Theorem (T} 

Proposition 4 Let 



UJi = Pl{z,t) +pi{z,t), UJ2 = P2{z,t) +P2{z,t), 



where pi and p2 are holomorphic functions of z which satisfy (14)- Let 
us consider the extended Moutard transformation of the operator Hq = dd 
defined by uJi, and let 9i he the image of uj2 under this transformation. Let 
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H = dd + U is obtained by the iteration of the Moutard transformation 
defined by 6i. Then 



U{z,z,t) = 2ddlogi\^{pip2 - P2P1) + J {{p'iP2 - PiP2)dz + {P1P2' - Pi'p2)dz) + 

, I I III III , ni I II II I \ , - - III - III - , n/ - II - I - I - ll\\ 

+ {P1P2-PIP2 + ^{PlP2 - P1P2) + PlP2 -Pl P2 + 2{pi P2 -P1P2 ))dtj. 

(15) 

Corollary 1 Let pi(z,t) and p2{z,t) be holomorphic functions in z which 



satisfy the equation (14)- 



Then the substitution of them into (15) gives a solution to the Novikov 



Veselov equation, which is rational in z, z, and t. 



Remark 4. It is clear that the analogs of (13) may be derived for 



all equations from the Novikov- Veselov hierarchy and hence the explicit 



solutions to them, in particular, of the form ( 13 ) may be constructed. 



4.2 (T-flows 



Let us denote by TIn the space of all harmonic polynomials p{z) + p{z) of 
the form 

p{z) 

and denote by Ti^j^ the subspace (Tq = 1. 



aoz^ + aiz^ H '^aN-iz + aN 



The flow (14) generates on TIn a linear flow 



af. = (N -k + 3){N -k + 2){N -k + l)ak-3, k = 0,...,N. (16) 

Since ctq is constant along the flow we restrict these equations onto Ti.^ and 
obtain a particular example of the linear system on cji, . . . , fiAr: 

& = A + Ba (17) 

which generates a dynamical system on the n-th symmetric product 5'"C 
of C. Indeed, ai, . . . ,an are the elementary symmetric polynomials in the 
roots zi, . . . , Zn of p{z): 

ai{zi, . . . , Zn) = -{Zi H \- Zn),...,(7n{zi,...,Zn) = {-!)"' Zi ... Zn 

and the integrable (even linear) evolution of a induces a dynamical system 
on S"'C. We call such a dynamical system on S"'C a cr-system. 
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In contrast to the Calogero-Moser systems, cj-systems do not describe 
the evolution of singularities of a solution to the NV equation (as the 
Calogero-Moser flow does for the KdV equation), i.e. the evolution of 



particle- type solutions [312] (see also ^2.2). In fact. 



the Calogero-Moser systems and a-systems are different in many re- 
spects. 



For example, let us consider the simplest flow induced by (14) on 7i-y. 
(Tl =0, £72 = 0, &3 = 6. 

The general solution is 

(Tl = ai, (T2 = a2, (73 = as + 6t 

However for a solution (xi(t), . . . ,x„(i)) to the Calogero-Moser system the 
k-th elementary symmetric polynomial ak{xi, . . . polynomial in t 

of degree k and especially any function Xk{t) is algebraic (see Corollary 3 on 
p. 118 in [2]). 

So we may present examples when the solutions to the Calogero-Moser 
system and to a a-system are related by a reparameterization of the time 
variable. For instance, such is the equilateral triangular solution to the 
Calogero-Moser system which corresponds to the solution 
-2^ [{x - xi{t)){x - X2{t)){x - X3{t))] of the KdV equation: 

Xk{t) = e'^^t, e = e^, A; = 1,2,3; 
and a solution to the cr-system corresponding to p{z, t) = + 6t: 

Zkit) = e^^/oi, e = e^, A; = 1,2,3. 
By Corollary [T] we have 



Given two solutions pi G and p2 € WjJ^ of (16), by a substitution 
of e^^^pi and e*'^2p2, where Ai and A2 are real- valued constants, into 
(15) we obtain solutions of the Novikov-Veselov equationj^ Therefore 
to each pair of solutions to (16) there corresponds an (S"^ x S'^)-family 



of solutions to the Novikov-Veselov equation. 

We think that due to their natural appearance and not only because 
of their relation to the Novikov-Veselov equation which we did establish, 
fj-systems deserve a special investigation. 



■^If we multiply, for instance, p\ by a constant G E then uj — pi + pi is multiplied by 
p and, since the Moutard transformation depends on uo taken up to a multiple, this does 
not change the transformation. 
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4.3 On the blow-up of solutions 



For simplicity let us write the formula ( 15 ) as 



U = 2dd\og^{pi,p2). 

(see the renomalization formula (11)). Due to the integration in ^ the 
function $(pi,P2) is defined up to a constant which we take real- valued 
to obtain a real-valued potential U{z,z). By Proposition |3] the potential 
V = V* equals 

It is easy to check the following: 

Corollary 2 The potential ^ (see Theorem^ is a stationary solution of 
the Novikov-Veselov equation. 

On the other hand one can obtain solutions of the Novikov-Veselov equa- 
tions which blow up in finite time: 

Theorem 4 The solution U{z,z,t) of the Novikov-Veselov equation ob- 
tained form the polynomials pi = i z'^ , p2 = z"^ i)z using PropositionU] 



and the integration constant C = —20 in (13) has the form 

with Hi = -12(^24^2 + Utx + 24ty^ + Uty + - 3x^y + 2x^ - 2x^y^ 
4x^y - 2x^y^ - QOx^ - 5xy^ - 4xy^ - 30x + y^ + 2y'^ - GOy^ - 30y ; 

H2 = {3x'^ + 4x^ + 6x2y2 + 3y4 _^ 4^3 + 30 _ mf. 

This solution decays at infinity as r^^ and obviously blows up for some t > 0. 

5 The cubic superposition formula 

In this Section we extend the algebraic superposition formula of L. Bian- 
chi |1] for three initial solutions oji, uj2, of the Moutard equation to the 
case of the Novikov-Veselov equation. In fact, already the extended Moutard 



transformation (13) should be considered as an "extended superposition 



formula" which produces a new solution 9 starting with two solutions a;, 99 



of the linear problem (12) with some initial potentials U, V, V*. In contrast 
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to the well-known algebraic superposition formula for two solutions of the 



Sine-Gordon equation, (13 1 is not algebraic and requires a quadrature. This 



is typical for (2 -|- l)-dimensional Backlund transformations. In this case an 
algebraic superposition formula exists for three initial solutions (cf. [9] for a 
detailed discussion of this phenomenon). 

Namely, as one may check, the following statement, obtained by L. Bi- 
anchi fl] for the case of the Moutard equation (pxy = u{x, y)ip, is also valid 
in our case: 



Theorem 5 If oJi, 0J2, CJ3 are three solution of (12) where U = Uq, V = 
Vq, V* = Vq is some given solution of the Novikov-Veselov equation 
and 61, O2 (transformed solutions of (12) with new Ui = Uq — 2dd(hiijOi), 
U2 = Uq — 2dd(hiuj2) and respectively transformed Vi, V* — another set of 
solutions of the Novikov-Veselov equation) are obtained from uj-^ via (13) 



then there exists a unique solution 9' of the 4th linear system (12) with the 
potentials U = U12, V = V\2, = V^i*2 such that 9' is connected to 9i, 92 



with the extended Moutard transformations (IS). This 9' is expressible with 
an algebraic formula 



= ^^(^2 



A 



UJlUJi 



-UJ2UJ2, 



(18) 



where lo'^, ^2 are obtained from uj\, UJ2 according to (13). 



It is useful to explain this statement and the "extended superposition for- 



mula" (13) by commutative rhombic and cubic diagrams ("Bianchi super- 



position cube" ) shown on Fig. [T| 

The rhombic diagram in the left part of Fig. [T] shows two solutions uji 
and UJ2 of ( 12 ) with the initial potential Uq and the respective Vq, Vq (which 
form a solution of the Novikov-Veselov equation ([T]) ) . Those tOi produce two 
new potentials Ui, U2 (and the respective Vi, V*) — two new solutions of the 



Novikov-Veselov equation (this is step one). The formula (13) then gives us 
(a set of) 9i — a solution of (12) with the potentials Ui, Vi, Vi as well as 
92 = —^01 — a solution of (12) with the potentials U2, V2, V2 . Using either 
9i or 92 we can find the potentials U12 = Ui + 29(9 log 0i = t/2 + 205 log 02, 
V12 = Fi +292 log 01 = ^2+252 log 02, V*2 = V^+2d'^\og9i = ^2* +2^2 log 02 
shown as the last fourth circle of the rhombic diagram (this is step two). As 
we see, this second step requires a quadrature and produces a free integration 
constant in the final potentials U12, V12, V*2. 

If we have done those two steps for each of the three pairs {uji,uj2}, 
{ljijLUs}, {102,(^3} shown on the right cubic diagram (and we have accumu- 
lated also three additional constants of integration), the Bianchi algebraic 
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Figure 1: Rhombic and cubic superposition formulas for (2 + l)-dimensional 
integrable systems 



superposition formula (18) gives a solution 6' of (12) for the set C/12, V12 



of the potentials obtained on step two. 



The formula ( 18 ) was proved in [3] for the case of the Moutard equation 
xy = U{x, y)uj. A straightforward computation shows that it preserves the 
dynamics of the potentials defined by ([T]), as Theorem [5] states. 
Using (12) and (18) one can see that 



• if we start with Uq = Vq = V* = and take arbitrary harmonic 
polynomials as the initial solutions Ui, i = 1, . . . ,N, then the poten- 
tials obtained on the Grst two steps are rational functions of {x,y,t). 
All subsequent potentials obtained on the next steps will be rational 
solutions of the Novikov-Veselov equation as well. 

• If we know some large set of solutions of the Schrddinger equation 
{dd + Uo)'tp = for the initial potential Uq, then using Theorem^ we 
can construct the corresponding solutions of {dd + Ui2)ip = using 
quadratures (on step two) if we simply set L03 = ^ in the superposition 



formulas (12) and {18). 



• This proves that we obtain potentials U12 integrable on zero energy 
level if we start from some integrable Uq, for example C/q = 0. 
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6 A remark on periodic and quasiperiodic solitons 



The similar construction can be applied to the case when the initial functions 
uji and UJ2 satisfy the Schrodinger equation with the constant potential u = 

(-A - k'^)uJi = (-A - k^)uj2 = 0. 

It appears that the second iteration may give an integrable potential which 
is periodic and smooth. Let us expose one such an example. 
Let 

oji = sin kx, UJ2 = sin(ax + by), + = k'^, 

and the Moutard transformation defined by oji maps UJ2 into 6i (modulo 
multiples of ^): 

b f cos(ax + by + kx) cos(ax + by — kx) ^ ^ 



2 sin kx \ a + k a — k 

with C = const and the new potential is 

_ , 9 cos^ kx 
u = k^ + 2 — ^— . 

sin kx 

Let us iterate the Moutard transformation using 9i as the generating solution 
and obtain a new potential 

^=-u + 2 ^^^^'+^^^^^ =k^-2A log(c^i0i) 

which is nonsingular for C large enough. An example of such a periodic 
potential u and the corresponding solution ipi = l/9i for the constants 
C = 3, a = 0, 6 = 1, A; = 1 is shown on Fig. 8, 9, page [22} 

As we see by iterating the Moutard transformation one may construct 
integrable two-dimensional non-singular periodic and quasiperiodic (when 
a/b is not rational) potentials. The complete basis of solutions for the 
constructed potentials u = —41112 may be obtained by applying (18| with 
LOs = exp (i(px -|- qy)), + q^ = k"^. The spectral properties and the NV 
evolution of such potentials deserve an additional study. 
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Figure 2: The potential ([7|. 



Figure 3: The solution '0i in Q. 




Figure 7: The solution ^2 in (10) 




